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Abstract 

The bb spectrum is calculated with the use of a relativistic Hamiltonian where the gluon-exchange 
potential between a quark and an antiquark is taken as in background perturbation theory. We 
observed that the splittings Ai = T(1D) — ^(lP) and other splittings between low-lying states are 
very sensitive to the QCD constant Ay(ri{) which occurs in the Vector scheme, and good agreement 
with the experimental data is obtained for Av(2-loop, nf = 5) = 325 rfc 10 MeV which corresponds 
to the conventional A M g(2-loop, n f = 5) = 238 ± 7 MeV, a s (2-loop, M z ) = 0.1189 ± 0.0005, 
and to a large freezing value of the background coupling: a cr it(2-loop, q 2 = 0) = a cr i t (2-loop, 
r — > oo) = 0.58 ± 0.02. If the asymptotic freedom behavior of the coupling is neglected and an 
effective freezing coupling a s t a tic = const is introduced, as in the Cornell potential, then precise 
agreement with Ai(exp) and A2(exp) can be reached for the rather large Coulomb constant astatic = 
0.43 ± 0.02. We predict a value for the mass M(2D) = 10451 ± 2 MeV. 

PACS numbers: ll.15.Tk, 12.38.Lg, 14.40.Gx 
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I. INTRODUCTION 



Recently the CLEO collaboration has discovered the first stable D-wave state in bot- 
tomonium with a mass M(1D) = 10162.2 ± 1.6 MeV which is consistent with the J = 2 
assignment [J] The observation of this state is also a great success of the potential model 
(PM) since just in the framework of the PM approach the strategy how to observe the ID 
state was developed and correct values for the radiative transitions and the M(1D) mass 
were obtained A study of the orbital and radial bb excitations which have relatively 

large sizes is also very important for a better understanding of the fundamental interac- 
tion in the infrared (IR) region. At present different conceptions about the value of the 
strong coupling in the IR region exist. In lattice QCD one finds for the coupling constant 
in the static potential, parametrized as a linear plus Coulomb potential, the small value 
^static = 0.23 (n/ = 0) — 0.30 (nf = 3) ^J. In phenomenological potentials the Coulomb 
constant is larger but spreads over a wide range from a sta tic ~ 0.33 — 0.39 0, 0] up to the 
rather large values 0.42 — 0.45 adopted in recent calculations 0, El- If the asymptotic free- 
dom (AF) behaviour is taken into account, then an even larger freezing value of the physical 
coupling, a s (l GeV) ~ 0.9 ±0.1, was determined from the analysis of the hadronic decays of 
the r-lepton 9j]. Our point of view discussed here is that after the discovery of the T(1D) 
the bottomonium spectrum already contains enough information to determine the freezing 
value of the coupling constant in an unambiguous way. 

Among the excited states in bottomonium the most precise information about the static 
interaction at large r can be extracted from an analysis of the IP, 2P, ID, and 2D states, 
because they have large sizes, lie below the open flavor threshold, and have no hadronic 
shifts. We shall show that instead of the absolute masses it is more convenient to use 
the splittings between orbital excitations since they are less sensitive to the choice of such 
parameters as the 6-quark mass m and in the mass differences the relativistic corrections 
(about < 20 MeV) are partially or completely cancelled. Note that for these states the spin- 
averaged masses M cog (nL), and therefore also the splittings between them, are known with 
great accuracy, ~ 1 MeV. For the nS states (since the states rj b (nS) are still unobserved) the 
accuracy of the spin-averaged masses M cog (nS) as well as the splittings 2S-1S and IP-IS is 
only about 10 MeV. 

From experiment we take two splittings between the spin-averaged masses for orbital 
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excitations: 



Ai = M cog (lD) - M cog (lP) = 262.1 ± 2.2 (exp)+J (th) MeV, 

A 2 = M cog (2P) - M cog (lP) = 360.1 ± 1.2 (exp) MeV, (1.1) 

which are measured, as well as the masses M cog (lP) = 9900.1 ± 0.6 MeV, M cog (2P) = 
10260.0 ± 0.6 MeV, and M(1D 2 ) = 10162.2 ± 1.6 MeV with high accuracy. However, for 
M cog (lD) one needs to take into account the very small difference between M cog (lD) and 
M(1D 2 ) coming from the fine structure splittings: M cog (lD) = M cxp (lD 2 ) + <5fs- in Ref. jjj 
5fs was found to be small lying in the range < 5fs < 1 MeV so that <5fs was included in 
the theoretical error in Eq. (Jl.lj) . 

In this paper we show that the splitting Aj weakly depends on the kinematics, the string 
tension a and the pole mass of the b quark, being at the same time very sensitive to the 
gluon-exchange potential used. Therefore Ai (as well as A 3 = M(2D) — M(2P)) can be 
used as a probe of the freezing value of the strong coupling in the IR region. 



II. HAMILTON! AN FOR SPINLESS QUARK AND ANTIQUARK 

In calculations of the bb spectrum we use the relativistic (string) Hamiltonian Hr which 
is derived from the Fo ck -Fe y nman- Schwinger (FFS) representation of the meson Green's 
function in QCD [l2|. This Hamiltonian and the calculated meson masses M(nL) 

contain only fundamental quantities: the current (pole) mass of the b quark, the string 
tension a, and the QCD constant Aj^g(nf) . It is known that in QCD, besides the current 
quark masses, only one "external" mass scale must occur, say the QCD constant Aj^g, while 
the string tension has to be expressed through this scale. Such a relation is established in 
Lattice QCD, however, in analytical approaches, including BPT, a relation between a and 
A is still not derived. 

The Hamiltonian is derived using several approximations. 

- First, the bb states below the open beauty threshold are described in the one-channel 
approximation, i.e., the creation of an additional quark-antiquark pair is neglected. 

- Second, the spin-dependent terms in the QQ interaction are considered as a perturbation. 

- Third, the minimal area law is used for the vacuum average of the Wilson loop. Then the 
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Hamiltonian has the following simple form: 



H R = UJ +^L + r- + V st (r) (2.1) 

00 UJ 



where p 2 = p 2 + L jr 2 and the "einbein" variable u appears when in the Green's function in 
the FFS representation one goes over from the proper time r to the actual time t: 2lu = $ L . 
According to the quantization procedure the operator u is determined from the extremum 
condition Q, 

dH R 



(Xc 0=> U =y/p* + m\. (2.2) 

Therefore the Hamiltonian (|2.1|) reduces to an expression which does not explicitly depend 
on u. 



H R = 2Jp* + m 2 + V st (r). (2.3) 



The kinetic term in Eq. (|2.3|) was postu 
the relativized potential model 



is postulated many years ago Q] and successfully used in 

0,0, HQ HQ, 



however, the derivation of this term 



directly from the QCD meson Green's function was done for the first time in Ref. Q|. The 
Hamiltonian (j2.1j) or ()2.3|) also implies a definite prescription for m q while in the PM the 
quark mass m q is usually considered as a fitting parameter. 

Let us consider two limiting cases. The first one occurs when the nonperturbative (NP) 
contribution dominates in the QQ potential (at relatively large distances) and the QQ 
gluon-exchange term for excited states can be considered as a perturbation [Lj]. Such 
an approximation is valid for light mesons {m q = 0). Then by derivation the mass m q 
in Hamiltonian ()2.1|) coincides with the Lagrangian current mass (the difference between 
current and pole masses can be neglected in the approximation considered): 

m q = m q (m q ), V(r) = V™ p (r) = or. (2.4) 

One must also take into account the NP self-energy correction to the quark masses [l7|. 
however, in bottomonium the NP self-energy contribution to H R appears to be compatible 
with zero (~ 3 MeV) and can be neglected Q]. For the linear potential the splitting Ai 
between the orbital excitations ID and IP turns out to be ~ 160 -r- 170 MeV, i.e., much 
smaller then the experimental value (jl.ljl . 

A different situation takes place at small distances where the perturbative QQ interaction 
dominates while in first approximation the string interaction can be neglected. Then due to 
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perturbative self-energy corrections the pole mass m q appears in the Hamiltonian ()2.1|) and 
for a consistent description m 9 (pole) has to be taken in the same n-loop approximation as 
the perturbative QQ potential: 

™, = <-(n- loop), yP (r) = _| °v(°-K>°P) , (2 . 5) 

For low-lying bb states this type of calculations was performed in Refs. (with the use 
also of the l/m q expansion of the kinetic term in Eq. ()2.3j) ). 

The perturbative static potential in Eq. (|2.5jl contains the vector coupling ay(r) in coor- 
dinate space. Its relation to the vector coupling cny(q) in momentum space and to the QCD 
strong coupling constant a s (q) in the MS scheme (up to three-loops) was studied in detail 
in Refs. [20] where the constant Av(Ar), which determines the perturbative vector coupling 
ay in momentum (coordinate) space, was expressed through the QCD constant A^: 

A v (nf) = A M g(n f ) exp (j^- j 

A R (n f ) = A v *, (2.6) 

where a± = ^ — ^§ri{, (3q = 11 — |n f , and 7^ is the Euler constant. 

If one takes the conventional value A^(2-loop) = 216±25 MeV, which corresponds to the 
"world average" a s (2-loop, M z ) = 0.1172 ± 0.0020 Q, then from the definitions Eq. (l2~oT) 
it follows that 

A(, 5) (2-loop) = 295 ±34 MeV, 

A^ 5) (2-loop) = 525 ±61 MeV. (2.7) 

Thus both Ay and A R appear to be essentially larger than A^ and therefore the pertur- 
bative description (where Apr C 1 or r < 0.3 fm), in a strict sense, can be applied only 
to the ground state T(1S). (The validity of the perturbative static potential in quenched 



approximation was discussed in Refs. |22l. |23j.) Since the orbital and radial excitations in 



bottomonium have rather large sizes, e.g., -R(IP) = 0.4 fm, R(1D) = 0.5 fm, R(2P) = 0.6 
fm, i?(2D) = 0.7 fm, for them the QQ gluon exchange cannot be described in the framework 
of PQCD alone and a NP modification of a s (ay) at large distances (IR freezing) should be 
taken into account. In our calculations we shall use Ay'* from the range (|2.7j) . 
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III. STATIC POTENTIAL 



To describe the bb spectrum as a whole we shall use here the static QQ potential as it is 
defined in background perturbation theory (BPT), where the influence of the background 
field on the gluon exchange is taken into account 
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2J: 



V st = ar + V B (r), (3.1) 

V»(r) = (3.2) 
3 r 

Here the background coupling a-p.(r) is defined through the Fourier transform of the potential 
Vb(qO in momentum space [23]; it gives 

«B(r) = - / — sin(gr)a B (g), (3.3) 
n Jo q 

where the background coupling in momentum space a^iq) is defined over the whole region 
< q < oo. For example, in two-loop approximation 

. 4vr 1 / /?ilnt B \ 

«B(g) = — — 1 1 - ^2 1 



A) *b \ Po h 
q 2 +_Ml 
A? 



t B = In ^2-^, (3.4) 



v 

where the background mass Mb under the logarithm is determined by the lowest hybrid 
excitation of the string. This mass cannot be considered as an additional (fitting) parameter 
since Mb itself can be calculated either on the lattice or with the use of the corresponding 
Tamiltonian for a hybrid 121. We take here the value of Mb which was obtained in Ref. 

~i . n 

23| from a fit to the static potential on the lattice at small distances |22j], 

M B = 1.00 ± 0.05 GeV. (3.5) 

Note that the most transparent way suggested to derive the expression ()3.4j) for a-Q^q) is to 
consider a large number of colors, iV c ^> 1 [ll] . 

The background coupling in momentum space a^iq) has an important feature-the correct 
PQCD limit at q 2 3> M| (as in Refs. J^) and therefore the constant Ay in the Vector- 
scheme ()3.4j) used in our calculations, is determined by the conventional Aj^g(rif) according 
to the relation (|2.6|) . 

One should keep in mind that the additive form of the static potential ()3.1|) is automati- 
cally obtained in the framework of BPT in the lowest approximation when the interference 



terms are neglected. This form is well supported by the lattice data |25( where the static po- 
tential is taken as the sum of a Coulomb-type term plus linear term at separations r > 0.2 
fm. Note that both potentials satisfy the Casimir scaling law with an accuracy of a few 
percent j^. 

In the next approximation of BPT the interference between perturbative and NP effects 
appears in the form of the background mass Mb (|3.4|) which ensures the IR freezing of ckb(?)- 
One can see from Eq. (|3.4j) that the only effect of Mb is to moderate the IR behavior of 
the perturbative potential, whereby the Landau ghost pole and IR renormalons disappear, 
while the short distance behavior, as well as the Casimir scaling property, stay intact 
Note that the logarithm in the coupling (|3.4|) formally coincides with that suggested in Refs. 



27| with Mg = Am 2 where a picture with the gluon acquiring an effective mass m g was 
suggested. However, since a physical gluon has no mass, the parameter (2m g ) 2 needs to be 
reinterpreted in the correct way. 

Our calculations of the bb spectrum are done for the number of flavors rif = 5 when 
the QCD constant A™ is well known from high energy processes, A^> = 216 ± 25 MeV, 

IVlo IVLo 

which corresponds to a s (M z ) = 0.117±0.002 20] and correspondingly in the Vector-scheme 

(5) 

A^> = 295 ±34 MeV. Later in our calculations we show that the splittings Ai and A 2 appear 
to be in agreement with experiment only for those values of Ay^ (2-loop) which are close to 
the upper limit in Eq. ()2.7|) . 

A$ } (2-loop) = 325 ± 10 MeV, (3.6) 

which corresponds to 

A^(2 - loop) = 238 ± 7 MeV, a s (M z ) = 0.1189 ± 0.0005. (3.7) 

From the definition ()3.4|) and for Ay^ ()3.6|) the freezing (critical) value of the background 
coupling a-B^q = 0) = a^ij- — > oo) is expected to be 

a crit (n { = 5) = 0.56 ± 0.01 (M B = 1.0 GeV), 

a C ritK = 5) = 0.59 ± 0.01 (M B = 0.95 GeV). (3.8) 

Notice that from Eq. (j3.3|) a crit turns out to be the same in momentum and coordinate 
space. 
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For a value of Ay close to the lower bound the freezing value a cr it Eq. ()2.7|) turns out to 
be essentially smaller, e.g., for A^ = 190 MeV (A^ = 260 MeV) the value a crit = 0.46 is 
obtained. For such a choice of Ay" 1 agreement with experiment for the bb spectrum cannot 
be reached. 

For A^, Eq. (Q, the pole mass of the b quark can be determined according to the 
standard procedure (as in PQCD) and for the conventional value rh^ffib) = 4.23 ±0.03 GeV 
one obtains in two-loop approximation 

m 6 (pole) = 4.82 ± 0.03 GeV. (3.9) 

Thus all fundamental quantities used in our calculations: m h (pole), A v and the string 
tension a ~ 0.18 ± 0.02 GeV 2 (taken from the slope of the Regge trajectories for mesons) 
are taken in correspondence to QCD; they are determined in the narrow ranges ()3.6)1 and 

(EH). 



IV. RELATIVISTIC CORRECTIONS 

Here for the bb spectrum we try to calculate the splittings Ai and A2 defined in Eq. (jl.lj) 
with high accuracy. First of all it is of interest to compare the spin-averaged masses M(nL) 
in the relativistic case, where M(nL) coincides with the eigenvalue (e.v.) of the spinless 
Salpeter equation (SSE), 

(2^p 2 + m 2 + V Bt (r)j ij){nL) = M(nL)ifj(nL), (4.1) 

with the e.v. M(nL) in the nonrelativistic (NR) limit which is widely used in bottomonium, 

M{nL) = 2m + E{nL). (4.2) 

Here the "excitation energy" E{nL) is the e.v. of the Schrodinger equation with the same 
static potential as in Eq. (|4.1j) : 

V„ = or - l55fe>. (4.3) 
3 r 

The mass m = mfe(pole) and «g(r) will be taken below in 2-loop approximation. 

In Section we also present the bb spectrum calculated with the Cornell potential which 
has been successfully used for many years ^|: 

V c (r) = or - + Co, (4.4) 



with a st atic = const. Since in a sta tic the (AF) behavior is neglected, this coupling can be 
considered as an "effective freezing" coupling in the gluon-exchange potential. Notice that 
a-s,{r) in BPT appears indeed to be close to a constant (the critical value) already at rather 
small QQ separations, r > 0.4 fm; e.g. for n f = 5 the ratio aB( r )/o ; crit is equal to 0.83, 0.88, 
0.92, and 0.95 respectively, for r equal 0.4 fm , 0.5 fm, 0.6 fm, and 0.7 fm. 

In Table |U the spin-averaged bb masses below the BB threshold, the splittings 

Ai = T(LD) - T(IP), A 2 = T(2P) -T(1P), A 3 = T(2D) - T(2P), (4.5) 

and the relativistic corrections <5r = M(nL) — M(nL) are presented for the static potential 
in BPT Eq. (14 .3J) with a typical set of parameters taken from the ranges ()3.6|) and ()3.9|) : 

m = 4.828 GeV, a = 0.178 GeV 2 , A^ = 330 MeV, M B = 1.0 GeV. (4.6) 

The accuracy of our numerical calculations is ±1 MeV for the e.v. of the SSE and ±0.5 
MeV for the Schrodinger equation. 

From Table H] one can see that the relativistic corrections <5r give a contribution (~ 
10 — 22 MeV) to the absolute values of the bb mass, nevertheless the splittings Ai and 
A 3 = M(2D)-M(2P) for the SSE and in the NR case appear to be equal (within ± 1 
MeV). This result is practically independent of the choice of such parameters as a and m. 
Therefore A x can be used as an important criterion to distinguish between different choices 
of the strong coupling. It can be shown that for the S-wave states <5r is larger and turns out 
to be more sensitive to the choice of m and o \ 

In the splitting A 2 = M(2P) — M(1P) the relativistic corrections are partly cancelled 
with A^ R being 7 MeV larger than A 2 for the SSE. Therefore it is better to use the SSE for 
the calculation of the highly excited bb states. (For the Cornell potential we find practically 
the same values for 5 R .) 



V. ANALYSIS OF THE SPECTRUM 

To demonstrate the sensitivity of the bb masses and the splittings to the choice of Ay f ^ 
present in the coupling a^ir) we give in Table HP the e.v.s MinV) of the SSE for three 
different values of Ay in two-loop approximation taken from the range (|3.6|) [a = 0.178 
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TABLE I: The spin-averaged bb masses M(nL) for the SSE (jHJ), M{nL) (jOJ in the NR case, 
the splittings Ai, A2, and A3 (|4.5|) . and the relativistic corrections <5r (all in MeV) for the static 
potential in BPT (j3.1|) with the parameters (|4.6|) . 





M(nL) 


M{nL) 




State 


SSE 


NR 


Sr 


IS 


9470 


9487 


-17 


2S 


10022 


10042 


-20 


3S 


10368 


10393 


-24 


IP 


9900 


9909 


-9 


2P 


10266 


10282 


-16 


3P 


10555 


10577 


-22 


ID 


10157 


10166 


-9 


2D 


10458 


10474 


-16 


Ai 


257 


257 





A 2 


367 


374 


-7 


A 3 


191 


191 






GeV 2 , M B = 0.95 GeV). We denote these three sets as 

I m = 4.819 GeV A^ 5) = 300 MeV, 
II m = 4.830 GeV A^ 5) = 320 MeV, 
III m = 4.836 GeV A^ 5) = 340 MeV, 

From Table ILD one can see that the splitting A x increases with growing Ay^ (or A^g.) and 

(5) 

reaches the experimental value at the value A v rj 325 ± 10 MeV which corresponds to 

A^(2 — loop) w 238 ± 7 MeV, a s (m z ) = 0.1189 ± 0.0005 (5.2) 

and gives rise to the critical value (M B = 0.95 GeV) 

«crit(n f = 5) w 0.58. (5.3) 

Thus, the conventional A^ (J5.2J) (close to the upper limit ()2.7j0 appears to be at the same 
time compatible with the large freezing value of the background coupling, giving rise to 
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Agi = 200 MeV, 

A i = 234 MeV > 

A® = 249 MeV. (5.1) 



TABLE II: The bb spin-averaged masses M(nL) (in MeV) for the SSE and the splittings Ai = 
M(1D) - M(1P), A 2 = M(2P) - M(1P), and A 3 = M(2D) - M(2P) for different values of m and 
(or A§L) given by Eq. ((BTTjl . The values a = 0.178 GeV 2 and M B = 0.95 GeV were kept fixed. 



State 


I 


II 


III 


experiment 


IS 


9484 


9473 


9468 


9460.3 ±0.3 (l 3 S"i) 


2S 


10023 


10023 


10023 


10023.3 ±0.3 (2 3 5i) 


3S 


10364 


10370 


10372 


10355.2 ±0.5 (3 3 5i) 


IP 


9900 


9900 


9900 


9900.1 ± 0.6 


2P 


10262 


10266 


10269 


10260.0 ±0.6 


ID 


10152 


10157 


10159 


10162.2 ± 1.6 (exp)ij;[} (th) 


2D 


10450 


10457 


10461 




Ai 


252 


257 


259 


262.1 ±2.2 (exp)lj^ (th) 


A 2 


362 


366 


369 


360.1 ± 1.2 


A 3 


188 


190 


192 





a good description of the bb spectrum and of the splittings A x and A 2 . Note that the 
critical value (J5.3)) is in striking agreement with a cr i t = 0.60 introduced in Ref. in a 
phenomenological way. 

This important statement does depend neither on the value chosen for a nor on the 
quark mass. To illustrate this fact we give in Tables IIHI and IIVI the splittings Ai and 
A 2 as well as the IP-IS and 2S-1S splittings for different values of a and two values of 
A v (n f = 5): A ( y } = 280 MeV which corresponds to A|| = 205 MeV and A^ 5) = 335 MeV 
which corresponds to = 245 MeV (all A's in two- loop approximation). 

Below we would like also to show the dependence of the splittings on the value of the string 
tension. We take a in the range 0.165-0.185 GeV 2 and fix the QCD constant A^ 5) = 280 
MeV. (This value of A^ 5) corresponds to A^ = 205 MeV and a{M z ) = 0.116 for M B = 1.0 
GeV.) Note that for such a value of Ay^ the freezing value, a crit = 0.488, appears to be 20% 
smaller than that in Eq. ()5.3|) . 

From the numbers presented in Table IIHI one can see that the splitting Ai is slightly 
increasing from the value 241 MeV up to 251 MeV while a is changing in a wide range, from 
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0.17 GeV 2 to 0.185 GeV 2 , being still 10 — 20 MeV smaller than the experimental number 
given in Eq. (jl.lj) . Also other splittings between low-lying states, e.g., M(2S)— M(1S) and 
M(1P)— M(1S), appear to be by 30 — 50 MeV smaller than the experimental numbers for 
any o from the range 0.17 — 0.185 GeV 2 . 

(5) 

Our calculations show that for a reasonable value of Ay one cannot reach agreement 
with experiment by variation of the parameter a only. Larger values of Ay^ (around 320 
MeV) are needed to get Ai « 260 MeV. 

TABLE III: The splittings Ai = M(1D)-M(1P), A 2 = M(2P)-M(1P), M(2S)-M(1S), and 
M(1P)-M(1S) (in MeV) for = 280 MeV (A® = 205 MeV) and different a (in GeV 2 ) 
(M B = 1.0 GeV, m b = 4.82 GeV). 



State 


a = 0.170 


a = 0.174 


a = 0.178 


a = 0.182 


a = 0.185 


experiment 


Ai 


241 


243 


246 


251 


251 


262.1 ± 2.2 


A 2 


347 


351 


356 


360 


363 


360.1 ± 1.2 


M(1P)-M(1S) 


396 


399 


402 


406 


408 


ss 430 a 


M(2S)-M(1S) 


514 


518 


523 


529 


533 


ps 550 a 



The situation appears to be different if a large value Ag^(2 - loop) = 245 MeV (a(M z ) = 
0.1194) is taken ( see Table HVJ). In this case the same three splittings turn out to be in very 
good agreement with the experimental numbers for a — 0.18 ± 0.02 GeV 2 . 



TABLE IV: The same splittings as in Table EH for A^> = 335 MeV (A§U2 - loop) = 245 MeV) 



State 


a = 0.170 


a = 0.174 


a = 0.178 


a = 0.182 


a = 0.185 


experiment 


Ai 


253 


255 


258 


261 


263 


262.1 ±2.2 


A 2 


360 


364 


368 


372 


376 


360.1 ± 1.2 


M(1P)-M(1S) 


423 


426 


430 


433 


436 


ps 430 a 


M(2S)-M(1S) 


544 


549 


554 


559 


563 


ss 550 a 



While the calculated 1D-1P, 2S-1S, and IP-IS splittings are in precise agreement with 
experiment for A^ 5) = 320 - 335 MeV (A§L = 234 - 245 MeV), the 2P-1P splitting in these 



12 



cases appears to be ~ 10 MeV larger than the experimental one. We expect that a small 
decrease of M(2P) (and also of M(3S)) can be reached taking into account the flattening of 
the confining potential due to additional pair creation. 

The sensitivity of the bb spectrum to the Coulomb constant in the Cornell potential is 
illustrated in Table IV! where three variants with fixed m = 4.80 GeV and a = 0.179 GeV 2 
are presented while astatic = constant is varied. From Table IVl one can see that only for the 
choice with large Coulomb constant a sta tic = 0.41 — 0.43, good agreement with experiment 
is obtained for all states with the exception of the ground state and partly of the 2S state. 

TABLE V: The bb spin-averaged masses M(nL) (NR case) for the Cornell potential (m = 4.80 
GeV, a = 0.179 GeV 2 ) with different values of 

^static and Cq (the masses and Co are given in MeV) 





a st = 0.3645 


age = 04071 


a st = 0.4263 






C = -74 


C = -32 


C = -16 


experiment 


M(lS)+5 AF (lS) a ) 


9421 


9434 


9451 


9460.3 + 0.3 (l 3 Si) 


Saf (IS) 


-34 


+37 


+79 




2S b *> 


10000 


9995 


9993 


10023.3 + 0.3 (2 3 Si) 


3S 


10337 


10344 


10347 


10355.2 + 0.5 ( 3 Si) 


IP 


9900 (fit) 


9900 (fit) 


9900(fit) 


9900.1 +0.6 


2P 


10249 


10259 


10263 


10260.0 ± 0.6 


ID 


10139 


10153 


10158 


10162.2 ± l-6±5;g 


2D 


10430 


10448 


10458 




Ai 


239 


253 


258 


262.1 ± 2.2(exp) ^-°(th) 


A 2 


347 


359 


363 


360.1 ± 1.2 


A 3 


183 


189 


192.5 





a ) The AF correction <5af(1S) is defined in Eq. Q5.4JI 
fe ) The AF correction is neglected for all states with the exception of the ground state. 



For the ground state T(1S) the AF behavior is very important both for the mass and 
for the wave function at the origin. Therefore one may take into account the difference 
between a sta tic m the Cornell potential (playing the role of the effective freezing coupling) 
and CfB,{r) considering it as a perturbation. The perturbative treatment gives rise to the 
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mass correction 

W13) = |( ° Mac ; aB(r) ) ■ (5.4) 

Owing to the AF effect this correction turns out to be positive (see second row in Table EJ 
for astatic > 0.39. However, a relatively large value for the AF correction, <5af ~ 40 — 80 MeV, 
for the ground state clearly indicates that the AF effect in the gluon-exchange potential is 
very important and cannot be treated accurately as a perturbation to the Coulomb potential. 

The meaning of the quantity a s t a tic m the Cornell potential can be understood if one 
introduces in BPT the effective coupling for a given nL state according to the following 
relation: 




TABLE VI: The effective coupling a e g{nL) (|5.5[) for the static potential with the parameters (|4.6|) 



state 


IS 2S 3S IP 2P ID 2D 


a e s(nL) 


0.39 0.41 0.43 0.47 0.47 0.50 0.50 



The values of a e s(nL) are given in Table IVTl from which one can see that for the low-lying 
states the effective coupling in BPT, a C R ~ 0.40 — 0.43 is indeed very close to the Coulomb 
constant used in the Cornell potential. 

From the spectra presented in Table IVl one can conclude that the splittings Ai and A2 
are in precise agreement with experiment for 

Static = 0.43 ± 0.02 (5.6) 

and this value appears to be 25% lower than the critical value for the background coupling 
states given in Eq. ()5.3|) and very close to the effective coupling given in Table EU 

For the best variants from Tables ITT1 and IVl the splitting A3 = M(2D) — M(2P) is practi- 
cally unchanged, A 3 = 191 ±2 MeV and therefore our prediction for the mass of the excited 
2D state is M cog (2D) « M(2D 2 ) = 10451 ± 2 MeV. 

VI. CONCLUSION 



From our analysis of the bottomonium spectrum we observe that the splittings between 
orbital excitations A 1; A 2 , and A 3 appear to be very sensitive to the freezing value of the 
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coupling in the gluon-exchange potential. The splitting Ai = M(1D)— M(1P) is of special 
importance since it is measured with the great accuracy 

For the Cornell potential, when the Coulomb constant a sta tic can be considered as an 
effective freezing coupling, precise agreement with experiment takes place for a sta tic = 0.43 ± 
0.02. 

For the QQ gluon-exchange potential with the strong coupling taken as in BPT, where 
the freezing value a cr it is fully determined by the value of Ay(rif = 5) in the Vector-scheme, 

(5) 

good agreement with experiment is obtained only for Ay = 325 ±10 MeV which corresponds 
to the conventional close to the upper limit: (2-loop)= 238±7 MeV and a s (Mz) = 
0.119 ± 0.001 and at the same time gives rise to a large freezing value, a cr i t ~ 0.58. 

We also predict that the mass of the as yet unobserved 2D state(s) is M(2Dj) « 
M cog (2D) = 10451 ± 2 MeV. 
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